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Abstract 

S" 

Q\ • The spin state of two magnetically inequivalent protons in contiguous àtoms of a molecule becomes 

ON , cntangclcd by the indirect spin-spin interaction (j-coupling). The degree of entanglement oscillates at the 

beat frequency resulting from the splitting of a degeneracy. This beating is manifest in NMR spectroscopy 
as an envelope of the transverse magnetization and should be visible in the free induction decay signal. 
The period (~ 1 sec) is long enough for interference between the linear dynamics and collapse of the wave- 
function induced by a Stern-Gerlach inhomogeneity to significantly alter the shape of that envelope. Various 
dynamical collapse theories can be distinguished by their observably different predictions with respect to 
this alteration. Adverse effects of detuning due to the Stern-Gerlach inhomogeneity can be reduced to an 
acceptable level by having a sufficiently thin sample or a strong rf field. 
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The measurement problem of quantum mechanics arises because the theory is silent as to when and 
how the transition takes place in dynamical processes between linear, deterministic Schròdinger evolution, 
and non-linear, stochastic wave-function collapse. One can sweep this problem under the rug only so long as 
experimental probes are unable to enter the transition region. Advancing techniques in mesoscopic physics 2 
Oh, make such an entry possible. The purpose of this paper is to suggest that we may already be able to 

investigate the transition region using well-developed techniques of nuclear magnètic resonance (NMR) . 

To see how we might go about distinguishing one theory of wave function collapse from another, let 
us suppose that we can arrange a dynamical transformation which, during a time í e , causes the following 
evolution of the spin state of two spin-1/2 partides: 



4)1 4> - 2-V3(| |}| j) + |T)lí))-l |)| |). (i) 

Suppose further that during this evolution an inhomogeneous (Stern-Gerlach) magnètic field is used as 
a measuring device, the result of the measurent being that the system collapses into one of the states i) | i) 
or | f)l í)- Spatial separation produced by a Stern-Gerlach field is proportional to the square of the time, 
and the spreading of the initially localized wave function is linear in the time, so there is a time t sg at which 
the ratio becomes unity, and we have a resolution of the two states. 

Consider now what will happen if we arrange to have: 

te ~ tsg • (2) 

Two types of entanglement oceur: the spin-spin entanglement oceurring in (1) and the entanglement of 
spin and space parts of the wave function as the latter separates into non-overlapping fragments under the 
infhience of the Stern-Gerlach field. One type of theory of wave-function collapse known as "spontaneous 
localization theory" (SLT) 3,4 attributes all collapse to spatial separation. In such theories the entangled 
state 2~ 1 / 2 (| |)| i) + | |)| j)) is not driven to collapse until the spatial separation has taken place. The 
situation is quite different, however, in another type of theory 5 known as "induced non-linearity theory" 
(INL). In such theories any linear hamiltonian which is able to distinguish the factorized constituents of an 
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cntangled state must induce a non-lincar term that brings about collapse. If INL theory were correct, the 
presence of the inhomogenous magnètic field would induce a non-linear term that would push the entanglcd 
spin-spin state towards collapse as that state is forming. Competition with the linear hamiltonian would 
thus occur all during the interval < t < t sg . Stopping process (1) at various times t in this interval leaves 
the system in diffcrent states for the two theories, and this difference can have observable consequences. 

We are now going to show that this scheme for testing a collapse theory experimentally can be implc- 
mented using two proton spins in contiguous àtoms A and A' of a molecule AA' . 

To begin with we must distinguish the sort of measurement we have in mind from conventional NMR 
spectroscopy 6 ' 7 . NMR spectroscopy probes the internal dynamics by detecting shfts in the Larmor frequency 
u in a strong static magnètic field B. If onc applies a transverse magnètic field of magnitudc b « B in a 
narrow frequency band including u, the dynamical effect is as if one moved to a frame rotating with the 
Larmor frequency. In this frame the static field appears much weaker, and the transverse field appears to be 
static. The internal dynamics then competes with the weakened magnètic interaction, and the eigenfrequen- 
cies of the resulting hamiltonian carry information about the dynamics. Since these freqüències also control 
the fourier distribution of the transverse magnetization, they can be detected by turning off the transverse 
field and fourier analyzing the free induction decay signal which appears as the spins relax. It is essential to 
maintain a high degree of homogeneity of the static field to prevent detuning of the Larmor frequency. But 
the testing of collapse dynamics requires that the field be inhomogeneous, and hence we cannot simply make 
use of small frequency shifts as a probe. 

To obtain such a probe we must find an experimental manifestation of cntanglcmcnt that is rclatively 
inscnsitive to detuning. Fortunately, as we shall see, such a manifestation will be found in a remarkable cor- 
relation between the degree of entanglement and a low frequency envelope of of the transverse magnetization 
(ETM). This correlation appears because both the entanglement and the ETM are produced by the same 
mechanism: a beating between two freqüències split by spin-spin interaction: 

The dynamical situation of the two protons that we are going to exploit is the following: A difference in 
their electrònic environments produces a different "chemical shift" of their magneto-gyric ratios and thereby 
of their Larmor freqüències. This effect is proportional to the B field strength and typically ~ lOppm (parts 
per million). There is a still smaller but measurable effect of the indirect spin-spin interaction (j-coupling) 
which is transmitted from onc spin to the other via a distortion in the electrònic wave function. The j- 
coupling, being a two-particle interaction, will entangle a spin-spin state provided that the state is not one 
of its eigenstates. Because of the chemical-shift the coupling to an external magnètic field is not simply 
proportional to a component of the angular momentum, and hence it does not commute with the j-coupling. 
Thus the state is not an eigenstate of the j-coupling and will become entangled. The frequency shift due to 
j-coupling is of order j w 1Hz (it is independent of the external field) and so the entanglement time t e will 
be of the order of seconds. 

Lct us compare this with the Stern-Gerlach separation time t sg . For dB/dz « lT/cm we have 

í- 1 - (3) 
~ h dz' {ò) 

where 7 is the average magneto-gyric ratio for the two protons and D is the molecular diameter. One finds 
that this is also of the order of seconds, so we have the required condition t sg w t e . 

Our next step is to examine the relationship between the entanglement and the transverse magnetization 

M=(^K/2|^) 2 + ^K/2|V) 2 ) 1 / 2 . (4). 

As shown in reference (5) there is a natural measure of entanglement: 

E(^)=2\detC\=2\C 11 C22-C 12 C2i\, (5) 

for the general state: 

\m) = Cu(í)| T)| T> + C 22 \ |)| I) + C 12 (í)| T)| 4) + C 21 (t)| 4)1 T>. (6) 
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One verifies that E ranges from for factorized states to 1 for maximally cntanglcd states. To compute E 
and M we must next determine the evolution of the state \ip{t)) with and without the INL contribution. 

The hamiltonian describing two protons with magneto-gyric ratios 71 ,72 in a magnètic ficld B is 
H B = (^<r (1) + ^<J (2) )-B = i(( 7l + 72 )<r+ (71 -72)^))- B, 

a = a^+a^\ Sa = <r« - a {2 \ (7) 

The term with 71 — 72 in Tíb is the chemical shift. We take the static contribution to be of magnitudc B 
in the z-direction and the rf contribution to be of magnitude b, rotating in the x-y plane with frequency lo. 
The j-coupling hamiltonian will be taken to be: 

H j =j*W·*M. (8) 

Going to a frame rotating with the frequency lo of the rf ficld, we obtain the effective spin part of the 
hamiltonian: 

H' = -{is<T z + d5a z + 2ja 1 ■ a 2 + \<r x }, 

v = ü> — u), d — ^(lüi — üj 2 ), \ = {b/B)ü). 

wi = jiB, 0J2 = 72-B, ü = ^(wi + W2). (9) 

Here we have assumed that b « B and so omitted the contribution of the chemical shift to the rf term. The 
parameter v, called the "detuning" , indicates the difference between the mean Larmor frequency of the two 
protons and that of the rf oscillation. Since we keep lo fixed, the variation in v comes from the difference in 
the value of ü for molècules in different parts of the samplc due to the inhomogeneity of the Stern-Gerlach 
field. 

The matrix 

h'= | x/0 J z,n 11a o7 I ( 10 ) 

represents H' in the basis: 

|1) = IT)IT), |2) = 11)11), |3) = |T)|I), |4) = ||)|T). (li) 
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The period of E can be deduced quite simply: For j '• — the odd terms in the characteristic polynomial 
for H' vanish, and so the eigenvalues have the form ±k , ±k\ with < k < K\. With j > these are 
perturbed, and in the parameter range above we find to first approximation: 

±k ±n Q + j, ±ki ±m - j. (12) 

Now observe that when we compute E we will encounter sums of products of two exponcntials of the form 
e lKt and e tK * where k, k' are any of the four eigenvalues given in (12). One sees that there is only one low 
frequency that can be obtained from a combination k + k', namely ±2j. In fact one obtains to lowest order 
in j: 

E = (l + is 2 /\ 2 y 1 \M23t)\=^te = ir/(2j), (13) 

provided that j « v. Thus one sees that the entanglement results from beating between the pair |k ± j\ 
or between the pair \k\ ±j|. If the initial state is | i)\ J.) one finds that the dominant contribution comes 
from the former. 
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These observations provide an important piece of information about the rangc of the detuning parameter 
v in which wc shall be intcrested. Because kq — > if v — > the beating is greatly reduccd for zero detuning. 
But the factor (1 + v 2 /A 2 ) -1 in (13) also reduces the entanglement as v becomes comparable to A, and so 
there is a middle range between j and A in which the entanglement is significant. 

We next perform a numerical computation to examine the relationship between M and E for a choicc 
of parameters that are expected to produce observable correlation. For protons at B = 1T (ü> « 40 MHz) 
a chemical shift of 10 ppm corresponds to d = 400 Hz. It is convenient to choose our encrgy or timc unit 
such that d = 1. We choose an rf field with b = 1G which gives A = 10 in our units, and we choose j = ÍHz 
which is 0.0025 in our units. We choose v = 5 as typical for the range 1 < v < X. 

In Figure 1 the transverse magnetization M is plotted in the absence of j-coupling so that there is no 
entanglement. In Figure 2 with j = 0.0025 the disentanglemcnt (1 - E) falls from unity to a minimum and 
then rises back to unity in approximately 1.5 sec, and the envelope of the transverse magnetization is visibly 
correlated with it. 



We next determine how the envelope is altered if, as the INL theory asserts, there is competition between 
the collapse mechanism and the linear dynamics during the interval < t < t e when t e w t sg . To do so we 
must first apply the INL theory of reference (5) to deduce the non-linear term appropriate to the experiment 
undcr consideration: The strength r\ of the term is dctcrmined by the reciprocal of the time t sg required 
for the measuring device to recognize the factorized constituents of the spin-spin entangled statc. Having 
arranged that t sg w t e , it follows from (13) that we must choose r\ w 2j. Since the Stern-Gerlach effect 
acts primarily on the states | |)| f ), | j}| J.) we may let the non-linear term act only in the subspace of the 
four dimensional spin space spanned by these two states. The stochastic element of the theory is introduced 
by a random fluctuation in the sign of the non-linear term. This also occurs in times ~ t sg and, since the 
proposed experiment is confined to < t < t sg , we may ignore this fluctuation. 

The non-linear, stochastic Schròdingcr equation of reference (5) as applied to the experiment under 
consideration can now be reduced to: 
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(15) 



The factor é 1 ar 9detc becomes indeterminate for det C — > 0, i.e. when the state factorizes. At this point 
the non-linear term turns off, which means that we may evaluate the expression as zero when it becomes 
indeterminate. The complex conjugation operation in the last term is an essential feature of the INL theory, 
the significance of which will be discussed below. 

The task of integrating (15) is greatly simplificd in our parameter range by the observation that arg det C 
is very nearly a constant during < t < t sg . The reason is the following: Suppose first that 77 = 0. If we 
compute C in the basis where H' is diagonal, we find that arg det C = tt/2 at all times. Since j is small, the 
transformation to that basis can be approximated by a tensor product of two one-particle transformations, 
and such transformations leave detC invariant. Thus for 77 = we conclude that arg det C remains close to 
7r/2 for < t < t e . Since 77 w 2j « 1 in our computations, we expect that this conclusion remains vàlid 
when the non-linear term is present. Suppose tentatively that this is the case. Then the four equations 
(15) together with the four complex conjugate equations are a linear system of eight equations which can be 
readily solved. We can then check the validity of our approximation by self-consisteney, i.e. computation of 
arg det C from the linearized solution and verification that it is indeed approximately tt/2. 

Figure 3 gives arg det C for 77 = (the near step function), the other curve being a plot of arg det C 
calculated in the linearized theory for r\ — 2j. Evidently the approximation is self-consistent in the first half 
of the interval, though deviations from the predictions of the linearized equation may be expected in the 
second half. 



4 



Figure 4 shows the transverse magnctization for r\ = 2j as calculatcd from the linearized equation. Its 
salient feature is a strong depression of the envelope as time approaches midpassage. This feature is not an 
artifact of the linear approximation because, as we have shown, that approximation is accurate during the 
first half of the period. 

Thus we arrive at a principal result of this paper: INL theory predicts a descent of the envelope of the 
transverse magnetization at times approaching the middle of the entanglement period. 

This descent will be found so long as the detuning is in the range 1 < v < A. If we assume that the 
gradient of the Stern-Gerlach field is dB/dz = 1 T/cm, we will have such a range over 90% of the samplc 
if the thickness is 1 micron. This restriction can bc rclievcd by raising the strength of the rf field. Thus 
if b = 10G the sample size can be increased to 10 microns with a corresponding increase in the effective 
signal. Ideally an experiment should first verify the presence of an envelope in the transverse magnctization 
that has the form of Figure 2 in a homogeneous B field, and then ascertain whether or not it deforms to 
that of Figure 4 when dB/dz is of the order lT/cm. The size of the deformation will place a bound on 7/, 
with r\ << 2j indicating failure of the INL theory. The absence of any change from Figure 2 when B is 
made inhomogenous would be consistent with spontaneous localization or with the commonly held view that 
collapse takes place instantaneously when one of the two factorized constituents is detected. 

To conclude we shall comment on the physics undcrlying the envelope distortion produced by the INL 
collapse mechanism: It was shown in reference 5 that the INL theory predicts a small CP violation in the 
decay of neutral K mesons which, with a very coarse model, comes within 20% of the experimental value. 
Neutral kaons are cntanglcd states of ds and sd quarks, and the non-linear term is induced by the semi- 
leptonic decay modes which can distinguish the factorized constituents. The source of the CP symmetry- 
breaking is the complex conjugation operation in the INL theory. Let us now show that the strong distortion 
of the transverse magnetization seen in Figure 4 is also a consequence of the complex conjugation operation 
in (15): 

Let the matrix appearing in the non-linear term be denoted T. It is anti-symmetric and has non- 
vanishing elements only in the (12) and (21) position. When we transform to the eigenbasis of H' it will 
be seen that T — > T' in which the latter is still anti-symmetric and has small off-diagonal clements except 
in the positions (12) and (21. In this basis the linear term produces oscillations e~ tat with a one of the 
four eigcnvalucs k — j, — «o — 3, K i + 3, an d — «i + j- Because of the form of T' just noted, the non-linear 
term produces a large coupling between the first two only. Were it not for the complex conjugation in the 
non-linear term, this would have little effect because of the large frequeney mismatch (2«o) bertween kq — j 
and — k — j. But, because of the complex conjugation, the coupling is between k — j and — (— k — j) so 
that the mismatch is only 2j. We thus have a near resonance between the linear and non-linear term and so 
obtain a large effect from a small perturbation that persists over the long entanglement period. 

The origin of the complex conjugation operation is the time-reversal operator, the oceurrence of which 
in the INL theory expresses the fact that while there is no arrow of time in microscopic processes, that 
symmetry is broken whenever a measurement is registered. Thus, if the INL theory is correct, it is this most 
bàsic attributc of the measurement process that manifests itself both in the CP asymmetry described in 
reference (5) and in the ETM depression descibed in this paper. 
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Figure Captions 

Figure 1: Transverse magnetization for no entanglement (j = 0) and a homogeneous B field. 

Figure 2: Correlation between transverse magnetization and entanglement (j = .0025). B - field homogeneous. 

Figure 3: Arg Det C for 77 = (the near step function) and for 77 = 2j. 

Figure 4: Transverse magnetization when there is entanglement (j = .0025) and INL collapse (i] = 2j.) 
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